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A phenomenological analysis of the distribution of Wilson loops in SU(2) Yang-Mills theory is 
presented in which Wilson loop distributions are described as the result of a diffusion process on 
the group manifold. It is shown that, in the absence of forces, diffusion implies Casimir scaling and, 
conversely, exact Casimir scaling implies free diffusion. Screening processes occur if diffusion takes 
place in a potential. The crucial distinction between screening of fundamental and adjoint loops is 
formulated as a symmetry property related to the center symmetry of the underlying gauge theory. 
The results are expressed in terms of an effective Wilson loop action and compared with various 
limits of SU(2) Yang-Mills theory. 



I. INTRODUCTION 

Wilson loops play an important role in gauge theo- 
ries. They describe closed gauge strings. Their expecta- 
tion value characterizes the phases of Yang-Mills theories. 
Extensive calculations of Wilson loops in lattice gauge 
theories have demonstrated confinement and shown the 
detailed behavior of the static quark potential. They 
have also revealed the existence of an intermediate dy- 
namical regime where simultaneously confinement pre- 
vails and so-called "Casimir scaling" occurs, in which 
the string tension in each representation is proportional 
to the quadratic Casimir operator 0. Asymptotically, 
for large loops, Casimir scaling must break down and be 
replaced by complete screening of adjoint charges and 
screening of half-integer charges to fundamental charges. 
In lattice calculations for SU(3) Yang-Mills theory, sig- 
nificant deviation from Casimir scaling has not yet been 
observed while in SU(2), indications of a transition to 
the screening regime have been obtained @, Q . 

Although these results of lattice gauge calculations 
have not yet produced a breakthrough in the understand- 
ing of the dynamics of Wilson loops, they have raised 
new questions. In particular, the origin of Casimir scal- 
ing turns out to be quite mysterious. Casimir scaling is 
expected theoretically for loops small enough to be calcu- 
lable in lowest order perturbation theory, or may occur in 
related approximation schemes such as the lowest order 
term of the cumulant expansion of the Wilson loop expec- 
tation value 5j. However a generally accepted explana- 
tion for the observed Casimir scaling in the confinement 
regime does not exist. Two-dimensional Yang-Mills theo- 
ries are known to exhibit Casimir scaling for loops of arbi- 
trary size and one may invoke "dimensional reduction" [lj 
to relate Yang-Mills theories in two and four dimensions. 
However, in two dimensions no screening occurs, and 
therefore the limit of asymptotically large loops is funda- 
mentally different from the four-dimensional case. Wil- 
son loops are complicated non-local, composite objects, 
which makes understanding their properties in terms of 
the fundamental degrees of freedom difficult. Hence, it 
is useful to take a phenomenological approach. Based 



on the observation that Wilson loop distributions calcu- 
lated in lattice QCD coincide, to a high degree of accu- 
racy, with distributions resulting from a diffusion pro- 
cess, we introduce a phenomenology based on diffusion 
on the group manifold. The resulting theoretical descrip- 
tion contains Casimir scaling as a natural limit and leads 
to a natural generalization that incorporates screening- 
induced deviations from this limit. This approach leads 
to the formulation of an effective action for Wilson loops. 



II. WILSON LOOPS 

The object of our investigations is the distribution of 
Wilson loops of fixed geometry, 



p{uj) = (0\S(lu-W) |0) 



(1) 



in SU(2) Yang-Mills theory. Wilson loops describe closed 
gauge strings along a prescribed curve C 



W = -trP exp{ig <f dx^A^x)} . 
2 Jc 



(2) 



In the confining phase, the vacuum expectation values of 
Wilson loops of sufficiently large area, A, obey the area 
law 



(0\W\0) 



-aA 



(3) 



where a denotes the string tension. The expectation 
value of a rectangular loop with side T much larger than 
side R is determined by the interaction energy of static 
quarks, V(R) 



(0\W\0) 



-TV(R) 



(4) 



As illustrated in Fig. ^ f° r small sizes, the Wilson loop 
distribution is peaked at values close to 1, whereas, for 
large sizes, it approaches the Haar measure 



A — > oo , p(u>) — > — yl 

7T 



(5) 
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III. DIFFUSION OF WILSON LOOPS 

Our description of the Wilson loop distribution starts 
with the observation that for a fixed gauge field configu- 
ration, either from lattice QCD or an ensemble of merons 
or regular gauge instantons, the value of a Wilson loop 
of a given shape and orientation undergoes a Brownian- 
like random walk as the size of the loop is progressively 
increased. 

On a lattice, for example, one may realize such a ran- 
dom walk by enlarging a sufficiently large rectangular 
loop in a fixed configuration by increasing the number 
of links on each side by one. To the extent that fluc- 
tuations produce confinement, the new Wilson loop can 
be thought of as a small, essentially random, step away 
from the previous loop. Similarly, for a fixed ensemble 
of merons or instantons, a small increase in the size pa- 
rameter of the loop, p, such as the radius of a circular 
loop, gives a new loop whose value can be expected to 
differ from the previous one by a small, nearly random 
step. Figure |21 illustrates how the values of differently 
oriented Wilson loops centered around a common point 
in a single gauge field configuration randomly evolve and 
become decorrelated as the area increases by the order 
of 1 fm 2 . (In the lattice calculation, to damp the ultra- 
violet fluctuations, we used APE smearing Q with 10 
smearing steps, each of which replaced each link by the 
original link plus 0.5 times the sum of staples connect- 
ing to that link.) Each such trajectory as a function of 
area constitutes a single random walk, and an ensemble 
of such random walk trajectories for an ensemble of con- 
figurations then yields the distribution of Wilson loops 
shown in Fig^for three areas. 

With this motivation, we now explore the assumption 
that each trajectory is exactly described by a Brown- 
ian random walk, where the time parameter, t, is an 
unspecified increasing function of the size parameter, p, 
of the loop. The distribution of an ensemble of such 
loops, p(oj,t), is expected to be described by a diffusion 



equation. For the case of SU(2) considered in this work, 
the diffusion of the (untraced) Wilson loop occurs on the 
group manifold S 3 . Due to gauge invariance, the value of 
the (traced) Wilson loop can be identified with the first 
polar angle $ on S 3 . The diffusion process is indepen- 
dent of the two other angles specifying a point on S 3 . 
When applied to functions which only depend on i9, the 
Laplace-Beltrami operator on S 3 , 



A.,3 = D»d u 



\ S 3 = u- Uf, (6) 



with the covariant derivative reduces to 



1 d . 2 . d 
■ sin 



sin 2 ■& dtf 9d 



(7) 



The diffusion equation on S 3 thus reads 



l-A s ,)G(0 1 t) = ^-6(d)8(t)- (8) 
at I sin v 



In terms of the eigenfunctions and eigenvalues of A53 
If sin(n + l)i? 



ip 7 . 



7r sin $ 
the spectral representation 
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n=l 



, -A53 ijj n = n(n + 2)tp n , (9) 



7rsin $ 



0(t) J2 n shintfe-t™ 2 - 1 )* (10) 



is obtained. 

Up to volume elements, we identify the Wilson loop 
distribution (JIJ with the solution of the diffusion equa- 
tion 



p(cosi?,i) =sintf G(d,t). 



(11) 



This identification does not specify the connection be- 
tween the size of the loop and the time t. However, 



FIG. 2: Values of Wilson loops in different planes centered at the origin as a function of the area of the loops (in fm 2 ) for a 
single configuration in SU(2) lattice gauge theory (left) and for a configuration of 500 merons (right). 



it predicts the shape of the distribution if the expecta- 
tion value of the Wilson loop and therefore the time t is 
known. 

As Fig. ^ and the left panel of Fig. demonstrate, the 
shape of the Wilson loop distribution follows closely the 
distribution of particles carrying out a Brownian motion 
on S 3 for a wide range of sizes and aspect ratios of the 
loops. In particular, the change in shape from the peaked 
distributions for small loops to the equilibrium distribu- 
tion, uniform on S 3 for large loops, 

2 

t — > oo, G —>■ — 

7T 

is correctly predicted, consistent with Eq. J5j|. 



IV. HIGHER REPRESENTATION WILSON 
LOOPS 

The expectation value of an observable 0(d) for the 
distribution generated by the diffusion process on S* 3 is 
given by 



2 />7T 1 — ' 

(O0&)) = - sintfdtfO (tf) Vnsinntfe" 



(n J -l)t 



Here, we consider the observables associated with the 
Wilson loops in the (2j + l)-dimensional representation 
of SU(2), which coincide with the eigenfunctions @ 



1 M 
— -tr exp{2itf I 



2j 

sin(2j + l)x) 



11 = 



r^W- (12) 



(2j + l) sintf V2(2j + 1) 
The j-dependence of the expectation values 

(Wj) = e -^ +1)t (13) 



is given in terms of the Casimir operator of SU(2) so that 
the ratios satisfy Casimir scaling 



(14) 



In lattice gauge theory, the validity of Casimir scaling 
has been demonstrated for SU(2) 1] and SU(3) |2| Yang- 
Mills theory. Approximate Casimir scaling has also been 
observed for ensembles of meron configurations 0. De- 
viations from Casimir scaling for sufficiently large loops 
have been observed in jjj, Q . 

In addition to the fact that Wilson loop distributions 
resulting from diffusion on S 3 imply Casimir scaling, the 
converse is also true. The Wilson loop operators in the 
different representations of SU(2) are given by the eigen- 
functions of A53. Therefore exact Casimir scaling fixes 
all moments of the Wilson loop distribution up to an 
overall constant, implying that the Wilson loop distribu- 
tion satisfies a diffusion equation with the diffusion time 
being a function of the size of the loop. This equivalence 
actually applies more generally to SU(JV). In the spec- 
tral representation corresponding to Eq. i|10fl , the "char- 
acter expansion" of G, the eigenfunctions corresponding 
to Eq. @ are the characters of the respective represen- 
tations and the eigenvalues are proportional to the cor- 
responding quadratic Casimir operator (cf. 



V. SCREENING 

In the Yang-Mills confining phase, two regimes in 
the Wilson loop dynamics can be identified. For small 
and moderate loop size, the Wilson loop distribution 
closely satisfies a diffusion equation and Casimir scaling 
prevails. However, asymptotically, for sufficiently large 
loop size, complete screening of adjoint (integer) charges 
and screening of the half-integer charges to fundamen- 
tal charges must occur. Screening of adjoint charges has 
been found in SU(2) lattice gauge theory 0, Q and for 
meron ensembles an indication of screening has been ob- 
served. The fact that the signal in the screening regime 



4 




FIG. 3: Approximate Casimir scaling for Wilson loops on a 32 4 lattice with lattice constant a = 0.13 fm. Left: Lattice Wilson 
loop distributions for loops of different aspect ratios with sizes given in lattice units. The solid lines are obtained from Eq. (11 (H 
with time t determined by the expectation values. Right: Logarithm of Wilson loop as function of the area (in lattice units) 
for j = 1/2, 1, 3/2, 2 representations. The straight lines are scaled by the corresponding ratio of the Casimir operators. 



is so weak suggests that the screening processes can be 
accounted for by a weak perturbation of the diffusion pro- 
cess. Hence, we extend our phenomenological description 
by adding a drift term to the diffusion equation. To this 
end we replace the Laplace-Beltrami operator, Eq. 
by 

A s3 = D^d^ A s3 = + (D^V) + (6^)6^ , 

which accounts for the presence of forces through the po- 
tential V. As above, due to gauge invariance, the Wilson 
loop variable can be identified with the first polar angle 
and the diffusion process must be independent of the two 
other angles on S 3 . We therefore require 

V = V(6) ■ 

The solution to the Fokker-Planck equation 

1 



has the spectral decomposition 



sin 2 7? 



5(ti)5(t) (16) 



t) = 6{t) J2 iWi(O) e"*"* e v ^ (17) 



with the eigenfunctions satisfying 



sin 2 d 



d-d 



■ 2 , dV ~ 
■sin 0—^ 



-A n ?/v, 



(18) 

The last factor in Eq. JT2J arises from the nonhermiticity 
of the differential operator (|15fl . The eigenvalue equation 
possesses the zero mode 



and therefore the asymptotic distribution in the diffusion 
process is no longer uniform on S 3 . Rather the equilib- 
rium distribution is determined by the Boltzmann factor. 

In order to discuss the properties of the screening po- 
tential V(i}), we note that the Wilson loop operators 
(|12|l obey the symmetry relation for reflection around 
the equator of S 3 (i9 — n/2) 



W i (7r/2- 1 9) = (-) 2 ^ j (7?- 7 r/2) 



(20) 



This symmetry property distinguishes Wilson loops in 
integer and half-integer representations and must be re- 
sponsible for their different properties. In order to avoid 
destroying this distinction, the dynamics of the Wilson 
loops must respect the reflection symmetry, so the po- 
tential V(#) has to be symmetric 



V{-k/2-$) = K(i9-tt/2) 



(21) 



Under the influence of V , the positive parity "free" states 
ip2n [cf. Eq. JJjJ] will mix among each other and so will 
the negative parity states ip2n+i- The symmetry of V 
prevents mixing of even and odd states. The expansion 
of the Wilson loop operators (|12fl in terms of the eigen- 
functions ip n reads 



^(*)=EC^(*) e 



VW/2 



, v = 0,l. 



i=0 



(22) 

In general, integer representation Wilson loop operators 
will contain a non-vanishing admixture from the zero 
mode 

This component will dominate at infinite times t [cf. 
Eq. i|17|) ]. and the expectation values of the Wilson loops 
in integer representations will become independent of t 



(19) 



(W„(<?)) 



(23) 
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The half-integer Wilson loops will also acquire a common 
asymptotic behavior, with the i-dependence determined 
by the lowest negative parity eigenvalue [cf. Eq. l(T%|) ] 



(W„+i(*)) 



w 2n-\ 



(24) 



Thus, reflection symmetry of V properly ensures asymp- 
totic screening of the integer Wilson loops to adjoint 
loops and screening of half-integer loops to the funda- 
mental loops. The transition from Casimir scaling to 
screening is controlled by the strength of V. 

We expect the asymptotics of the Wilson loops in in- 
teger representations to depend on the linear size of the 
system. For large loops, screening requires the formation 
of two colorless states, so called gluelumps (cf. @). In 
turn, the interaction energy [cf. Eq. (@J] of static adjoint 
charges must satisfy a perimeter law, with the slope de- 
termined by twice the mass of these objects m g i. This 
change in the asymptotic behavior requires a time depen- 
dent potential 



V = V($,t) 



(25) 



which asymptotically decreases with time. For large t, 
from the area law for half-integer loops 

t~p\ 

The correct asymptotics for integer loops is reached, pro- 
vided 

fi~2m gl . (26) 

When dynamical fermions are present also half-integer 
Wilson loops will satisfy a perimeter rather than an area 
law. In this framework, this essential change due to 
fermions can be incorporated by admitting potentials V 
which violate the symmetry relation, Eq. (|21|l . In this 
case, one expects all the Wilson loops to approach a com- 
mon asymptotic behavior determined by the mass of the 
lightest meson. For heavy quar ks and circular loops the 
potential can be calculated |13|. with the result 



/' 1 



2m c , 



(27) 



with the quark mass m q . The quark induced potential V q 
is antisymmetric around the equator of S 3 and therefore 
turning on Vq mixes even states (ip2 n ) and odd states 



VI. WILSON LOOP EFFECTIVE ACTION 



(2) 



Our results suggest treating the Wilson loop variables 



W = cos a 



(28) 



as quantum mechanical degrees of freedom with the fol- 
lowing effective action 



S c fi[a] = 



dt{la 2 + U(a)}. 



(29) 



For a quantum mechanical system the path integral 

[•a(t)=ti 

g(t',t;ti',#) = / d[a] e - s M, 

Ja(t')=i?' 



d[a] = TT / sin 2 aj doj 

4 JO 



(30) 



satisfies the diffusion equation on S 3 



d_ 

dt 



A S 3 + lM g(t',t;&,#) 



5(0-0') 5(t-t>) 

sin v 



(31) 

In terms of g, the Wilson loop distribution is expressed 



p(tM) 



g(0, t; 0, f) J dflpo sin 2 ggg g(t, oo; fl, tf^) 
J di?oo sin 2 g(0, oo; 0, i?oo) 
Xo(<?) 



= ff(O,i;O,0) 



Xo(0) 



(32) 



In the evaluation of the integral we have used the spectral 
representation of the Green function g. At infinite times, 
only the ground state of the system 



(-As3 + ET)Xo(0)=O 



(33) 



contributes. It is straightforward to verify that the Wil- 
son loop distribution satisfies the Fokker-Planck equation 



|-A S3 ) p(0,i) = ^i-W(i), (34) 
at J sin v 



with the potential V appearing in — A53 [cf. (|I5[l ] given 
by 



(35) 



The potential V , which may in principle be obtained 
from a fit to the Wilson loop distribution, determines 
the potential U appearing in the Wilson loop effective 
action 



U(a) 



dV 
da 



cot a 



1 dV_ 
4 ~da~ 



1 d 2 V 



2 da 2 



(36) 



Reflection symmetry of V around a — ir/2 implies that 
of U and vice versa, provided the ground state wavefunc- 
tion xo is also even under reflections. The results of our 
investigations can be described by the following effective 
action of the Wilson loop variables depending on the size 
parameter p 

dpr-\p) |I( T(/a )0) 2 + C/ (a)| . (37) 



S e s[a] z 
The function 



(38) 
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remains unspecified. Irrespective of the choice of r, 
the Wilson loop distributions generated by this action 
are solutions of the Fokker-Planck equation (|34|) . In 
the limit U = 0, Wilson loops exhibit Casimir scal- 
ing. In this limit, the action possesses a shift symmetry 
a(p) — > a{p) + 6a, corresponding to translational sym- 
metry on S 3 . Clearly, any a-dependent V destroys this 
symmetry and thereby gives rise to violations of Casimir 
scaling. 

The effective action separates the dynamics of diffu- 
sion, with the associated phenomena of Casimir scaling 
and screening, from the Wilson loop dynamics governing 
the size dependence, with the appearance of perimeter or 
area laws incorporated in r(p). 



VII. CONNECTION TO YANG-MILLS 
THEORY 

In the following discussion, we will indicate the con- 
nection of our phcnomcnological findings to SU(2) Yang- 
Mills theory. We will restrict our considerations to cir- 
cular loops (radius p) and discuss first the limit of small 
loops. 

In the short-time limit of vanishing Wilson loop size, a 
regular potential U becomes irrelevant. The Wilson loop 
distribution, Eq. ljTU|) . 



p(0,t)»-r 
TTV 



dnnafan*e-"* = 5-^5 e-*"/« 

(39) 

is obtained by diffusion in _ZR 3 , the tangent space of 
SU(2). For small sizes, the Wilson loops can be calcu- 
lated perturbatively. To leading order, the SU(2) Wil- 
son loops coincide with the free U(l) 3 gauge theory. In 
U(I) gauge theory the Wilson loop distribution is given 
by(i? G M) 



Pu(i)(tf,P) 



with 



| j d[A] 6(0 -g £dx,A^)e 



S[A] 



exp 



Ann J 



9_ 

2tt 



f dx^ f dy 1/ D /J , JJ (x - y) . 
Jc Jc 



(40) 



(41) 



ic Jc 

The singularity in the Feynman propagator, 

1 



Dy.v{x) 



Att 2 x 2 



6uu 



can be regularized by Gaussian smearing of the point 
charge over distances ~ e perpendicular to the loop plane. 
To leading order in e 



(42) 



In SU(2) Yang-Mills theory, for small loop size, the Wil- 
son loop distribution, given by 



# 2 d-d = l[p X j {1) (<& i , P )dti i , 



(43) 



i=l 



is the distribution generated in a diffusion process in M 3 . 
From comparison with (|39|l . we conclude 



0, i = 



4 v / 2, 



P, 



(44) 



and therefore, for small Wilson loops the effective action 
is given by 



5 e ff[a] = 



V2e 



(45) 



In the limit of large loop sizes, the fundamental loops 
obey an area law which, in the Casimir scaling regime 
[Eq. iJTUJl] requires the identification 



t = 



(46) 



Thus for large loop sizes 

Assuming the logarithm of the Wilson loops to be a sum 
of constant, perimeter and area terms, the function r(p) 
respecting the two limits Eq. Ij44(l . Eq. (|46|l is given by 



-(P) = 



27TCT 



n -1 



4V2t 



(48) 



Finally, we compare the effective action with the Yang- 
Mills action. To this end, we choose the gauge in which 
circular loops appear as elementary degrees of freedom, 
and assume the loops are located in the (x, y)-plane and 
centered around the origin. We choose cylindrical coor- 
dinates p,tp,z,t and transform to the azimuthal gauge 
[l3| with the (^-component of the gauge field 



1 T 3 

A v (x) = a(p,z,t) — 

girp 2 



(49) 



being diagonal and independent of ip. With this defini- 
tion, the gauge field a is the Wilson loop [cf. Eqs. J2J, 
(|28p] . In field strength components F 3 a containing 
derivatives of A v , no commutator terms are present, e.g. 



F 3 pv = l p (^d p a-d v A 3 p 



(50) 



4V2vr e 



The effective action of the Wilson loop variables is de- 
fined by 

e" s ™[ Q l = / d[A',a']e- s ^ A '- a ^Y[6{a'{p,0,0)-a(p)^ 

(51) 



7 



In the Yang-Mills action 



Sym = / P dp dip dz dt 



1 



\2-K 2 g 



(d p a) 2 + (8 z a) 2 + {d t af 



C 



(52) 

'}• 



due to the ip independence of a, no bilinear mixing of 
the Wilson loop variables with the other fields occurs. 
£' contains the interaction of the Wilson loop variables 
with the other fields. The integration over the Wil- 
son loop variables in Eq. (|51l) includes, as in Eq. (|3U[) . 
the Haar measure. For comparison with the results in 
Eq. 1|45(1 and Eq. I|47|l we observe that the effective ac- 
tion of the Maxwell theory is obtained by dropping the 
Haar measure in Eq. (|5ip. integrating a' over the real 
axis, and disregarding the a'-independent term CJ in 
Eq. 15211 . Integrating out the Wilson loop variables at 
(z,i) 7^ (0,0) apparently modifies the volume element 
in the effective action. Propagation of the Wilson loop 
variables in Maxwell theory occurs in (3+l)-dimensional 
space-time and therefore actively involves the degrees of 
freedom located at {z,t) ^ (0,0). As Eq. (|4*7|) suggests, 
this is not the case for large Wilson loops in Yang-Mills 
theory. Here the integration over the degrees of free- 
dom other than a(p, 0, 0) renormalizes the bare action 
but does not change the volume element. The nearest 
neighbor interaction seems to be operative only for small 
loops. This can be interpreted as a consequence of the 
formation of a flux tube, which turns the effective action 
into that of a (l+l)-dimensional system. Indeed, evalua- 
tion of k, Eq. (|41|l . using the two-dimensional propagator 
D^vix) ~ (5 p „lnx 2 yields the effective action, Eq. l|4*7|) . 

Our last topic concerns the evaluation of the Wilson 
loop distribution [Eq. JTJ] in SU(2) Yang-Mills lattice 
gauge theory. In the strong coupling expansion the dis- 
tribution of the values of a square Wilson loop of area A 
is given by 



2 

p(cosz9,t) = — sinz9 

7T 



l + 4costfe- tln9 + (53) 
+ (4cos 2 0-l)e- 16 <^- 1 > In <' a + . 



where, in terms of the lattice constant 



t = A/a 2 



(54) 



We have kept only the leading contributions in 1/g 2 to 
the area and perimeter dependent terms. To this order, 
the distribution function gives rise to non-vanishing ex- 
pectation values for loops in the fundamental and the 
adjoint representation. The area term arises from tiling 
of the minimal surface enclosed by the loop, the perime- 
ter term from tiling the minimal surface of a rectangu- 
lar "tube" of volume ~ 4^/Aa 2 along the loop. As ex- 
pected, in the distribution the area and perimeter con- 
tributions are distinguished by the reflection symmetry 
in u> = cosd. Symmetric and antisymmetric compo- 
nents p(u), t) ± p(—oj, t) contain only terms which satisfy 



a perimeter and an area law respectively. For large loops 
therefore, p(u>,t) becomes an even function in lo. We 
also observe that, to leading order, strong coupling im- 
plies screening of Wilson loops in integer representations. 
Thus in strongly coupled lattice gauge theory, Casimir 
scaling is not a natural limit. Since Casimir scaling im- 
plies an area law for loops in every representation, tiling 
by including contributions from links transverse to the 
loop must be suppressed by dynamics beyond the strong 
coupling limit. We thus reach in lattice gauge theory the 
same conclusion as in the continuum theory discussed 
above. In the limit of Casimir scaling, the effective the- 
ory has to be two-dimensional. Following the line of ar- 
guments given above, we can reach this limit formally 
by assuming that links transverse to the loop (in the 1-2 
plane) are not affected by the presence of the loop. The 
distribution function [Eq. QJ] is thereby reduced to a 
functional integral over the links in the plane of the loop 
U\2 and planes parallel to it 



p(cos$, t) 



J d[U 12 ] e' s ^ u ^ S (cost? - itrlL ^2) 



fd[U 12 ]e- s ^ u ^ 



(55) 



To simplify the action 

Si2[Ux 2 ] = 5-tr P 12 (ni,n 2 ,n 3 ,n4) , (56) 

q z * — ' 

Iii,n2,n3,n4 

we assume that the loop affects the plaquettes P\ 2 only 
in 1-2 planes within a distance A from the plane of the 
loop and that the contributions to the action from all of 
these planes are the same. With these approximations 
one may construct a crude but suggestive model of the 
flux tube picture of confinement. With 

Sxi[Uxi\ w -^r^tr P Mni,n 2 ), (57) 



9{ 2) a 



ni ,n 2 



we arrive at two-dimensional Yang-Mills theory with cou- 
pling constant 



9(2) 



9 

ttA 2 ' 



(58) 



In two-dimensional Yang- Mills theory, the partition func- 
tion and various observables can be evaluated in closed 
form, cf. [lfl lll| and some recent applications thereof 
[l2^ . We obtain the following expression for the Wilson 
loop distribution 



p(cos$, t) 



o 00 



sin nd 



In (2/9) V 



(59) 



in terms of the modified Bessel functions depending 
on the two-dimensional (inverse) coupling constant 



2ir\ 2 

2 2~ 

9 a 



(60) 
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and t as given in Eq. (|54|l . In the continuum limit 



oo, 



h(2(3) 



4/3 



we finally obtain 



2 

p(cos i) w — n --iii 

^ n=l 



(61) 



(62) 



with the string constant in the fundamental representa- 
tion 



3g 2 
8ttA 2 



(63) 



In this way, we reproduce the distribution obtained from 
diffusion [cf. Eq. I|lt)(l ] and thereby confirm the intimate 
connection between Casimir scaling and two-dimensional 
dynamics. Whereas confinement in two dimensions is of 
kinematical origin, Casimir scaling is obtained only in the 
continuum limit. In the strong coupling limit (f3 — > 0), 
the n 2 — 1 dependence on the degree of the representation 
in Eq. (|fJ2|) is replaced by a linear n dependence. 



VIII. CONCLUSION 

In this work we have presented a phenomenological 
analysis of the distribution of Wilson loops in SU(2) 
Yang-Mills theory. We have shown that Wilson loop dis- 
tributions calculated in lattice gauge theories or obtained 
in model studies can be succinctly represented by distri- 
butions resulting from diffusion of Wilson loops on the 
group manifold. This description is applicable to small 
loops in the regime of asymptotic freedom as well as to 
large loops in the confinement regime. As preliminary 
results [l3| indicate, it applies equally well in the decon- 
fincd phase. 

Casimir scaling of Wilson loops in different represen- 
tations is a natural limit in this description. It is an 
exact property of Wilson loops if and only if the diffu- 
sion occurs in the absence of external forces. It reflects 
the translational symmetry on the group manifold. This 
equivalence of Casimir scaling and free diffusion remains 
true for SU(iV). 

Screening processes can be accounted for by a poten- 
tial term in the diffusion equation, or equivalently, a drift 
term in the Fokker-Planck equation. The crucial distinc- 
tion in the screening process between fundamental and 
adjoint loops is formulated in this framework as a symme- 
try property of the potential term. A reflection symmetry 
protects the area law of large fundamental loops from be- 
ing violated by screening processes. The initial condition 
in the diffusion process breaks this symmetry by singling 
out the point on the group manifold where the source 



is located. Consequently, the solution of the diffusion 
equation makes this symmetry manifest only at asymp- 
totic times, that is, for asymptotically large loops. This 
analysis suggests the existence of a similar symmetry to 
distinguish fundamental and adjoint loops in the underly- 
ing gauge theory. The associated symmetry transforma- 
tions must include a Z2 center reflection of Wilson loops. 
Unlike the Z2 symmetry related to Polyakov loops, this 
center symmetry exists irrespective of a compact direc- 
tion in space-time. 

The formulation of our results in terms of an effective 
Wilson loop action separates confinement and screening 
dynamics. Confinement is primarily connected to the 
properties of the kinetic term in the action. The existence 
of an area law is related to the effective two-dimensional 
structure of the kinetic term for large loop sizes. This 
remains true as long as the potential respects a reflection 
symmetry. 

The strength of the potential term controls the devi- 
ations from Casimir scaling. For a non-vanishing poten- 
tial term, Wilson loops in half-integer representations ap- 
proach a common asymptotic behavior and so do Wilson 
loops in integer representations. If the potential violates 
the reflection symmetry, as induced for example by dy- 
namical fermions, an area law can only show up, if at all, 
at intermediate sizes. Asymptotically, the distinction be- 
tween half-integer and integer representations in the size 
dependence of the loops disappears. 

Extension of this work to SU(3) and more generally to 
SU(iV) is of interest. In SU(iV), the starting point will 
be a diffusion equation for N — 1 "radial" variables with 
associated Jacobians (cf. These variables cor- 

respond to the N — 1 Wilson loops that are elements of 
the Cartan subalgebra. The symmetry of the potential 
is crucial for the distinction of confined and nonconfincd 
variables. For SU(3), in a different context, the choice of 
variables and the symmetry analysis has been carried out 
15]. It will be of interest to see whether the suppression 
of screening processes in SU(3) as compared to SU(2) 
can be related to differences in the center symmetry. For 
N > 4, a new element in the dynamics appears. The 
string tension of strings of different Zjv charge may be 
chosen independently of each other. As a consequence, 
Casimir scaling might be violated or new scaling laws 
might emerge as string or supersymmetric theories sug- 
gest and lattice calculations indicate 0. A phe- 
nomenological symmetry analysis of the potential term 
in the corresponding Fokker-Planck equation may also 
offer a new perspective on the dynamics of Wilson loops 
and help to clarify their large- N limit. 
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